Defect formation in inhomogeneous 2-nd order phase transition: theory and 

experiment. 



G.E. Volovik 

Low Temperature Laboratory, Helsinki University of Technology, Box 2200, FIN-02015 HUT, Espoo, Finland 

and 

Landau Institute for Theoretical Physics, Moscow, Russia 
(February 1, 2008) 



On 
On 
On 



The status quo in our understanding of defect formation during a rapid transition into the broken 
symmetry state in condensed matter and cosmology is discussed. An observation of vortex nucleation 
in neutron absorption experiments in superfluid ^He-B is interpreted in terms of defect formation 
during inhomogeneous coohng through Tc. Due to the temperature gradient in the locally heated 
region the superfluid phase transition occurs as a propagating front. The theoretical considerations of 
vortex formation at the propagating front are based on work by Kibble- Volovik, Kopnin-Thuneberg, 
and Aranson-Kopnin-Vinokur (AKV). 
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I. INTRODUCTION 

To produce a new vortex line in the vortex-free state 
of superfluid liquid is not an easy job. If the container 
is devoid of the remnant vorticity, which can be pinned 
by rough surface, the vortices are created only when a 
treshold Vc for the hydrodynamic instability of the super- 
flow is reached [Q. The thermal activation or quantum 
tunneling can assist the nucleation only in the narrow 
vicinity of the instability treshold, where the external 
perturbations, however, are more effective. In super- 
fluid ^He-B, because of the large size Vc of the vortex 
core, the region near the treshold, where thermal activa- 
tion or quantum tunneling can be important, is partic- 
ularly small, Vc — Vs ^ 10~^Vc' In a typical cylindrical 
container with radius R = 2.5 mm and height L = 7 
mm, rotating with angular velocity Q = 3 rad/s, the 
vortex-free state stores a huge amount of kinetic energy 
(1/2) / dV Ps{vs — ^ 10 GeV. This energy cannot 
be released, since the intrinsic half period of the decay of 
this metastable state is essentially larger than the proton 
life time. 

That cosmic rays can assist in releasing this energy by 
producing vortex rings, I first heard from my supervi- 
sor, professor lordanskii, in 1972. The natural scenario 
for that was thought as depicted in Fig.|l|. The ener- 
getic particle heats a region above the superfluid tran- 
sition temperature Tc. During the cooling the normal 
liquid in this region can continuously evolve to form the 
core of a vortex loop, which starts growing if the radius 
Ri) of the heated region is larger than the radius of the 
ring sustained by counterflow, i.e. \i \vs — Vn\ > Vd = 
{f<i/4:7TRi})\n{Ri)/rc)^ where hi is the quantum of circula- 
tion around the vortex. 

If the counterflow essentially exceeds this treshold, the 
evolution, which is most favourable for vortex produc- 
tion, leads to the closely packed vortex rings of the crit- 
ical size, which can further develop. This gives the fol- 
lowing estimation for maximal number of vortex loops. 




FIG. 1. Evolution of the overheated normal fluid region 
into the core of the vortex ring in the presence of the coun- 
terflow. View from 70's. 



which can grow further: N ~ {\vs — Vn\/vci)^. 

Experiments with irradiated superfluid ^He were 
started in 1992 in Stanford, where it was found that the 
irradiation assists the transition of supercooled ^He-A 
to ^He-B. In 1994 the neutron irradiation of ^He-B was 
found to produce a shower of quasiparticles in Lancaster 
and vortices in rotating ^He-B in Helsinki ||]. Energy 
deficit found in low-T Grenoble experiments indicated 
possib le formation of vortices in ^He-B even without 
rotation In Helsinki the observed number of vor- 
tices produced per one event showed both the treshold 
behavior and the cubic dependence at large rotation ve- 
locity: Above the treshold it was well approximated by 
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N ^ {\vs — Vn\/vci)^ — 1. This indicated that the nature 
has chosen some scenario, which produces the maximal 
possible number of vortices. What is the reason for that? 

The decay products from the neutron absorption reac- 
tion generate ionization tracks, the details of which are 
not well known in liquid ^He. At the moment we have 
two working scenaria of thermalization of the energetic 
particles: 

(i) The mean free path is long and increases with de- 
creasing of the energy. This can lead to a "Baked Alaska" 
effect, as has been described by Leggett A thin shell 
of the radiated high energy particles expands with the 
Fermi velocity vp, leaving behind a region at reduced T. 
In this region, which is isolated from the outside world by 
a warmer shell, a new phase can be formed. Such Baked 
Alaska mechanism for generation of new phase has also 
been discussed in high energy physics, where it describes 
the result of a hadron-hadron collision. In this relativis- 
tic case the thin shell of energetic particles expands with 
the speed of light. In the region separated from the ex- 
terior vacuum by the hot shell a false vacuum with a 
chiral condensate can be formed [^. This scenario pro- 
vides possible explanation of formation of the B-phase in 
the supercooled A-phase [|[. 

(ii) During thermalization the mean free path is less 
than the dimension of the region where the energy is de- 
posited and the temperature is well determined during 
the phase transition through Tc. In this case there is no 
Baked- Alaska effect: no hot shell separating the interior 
region from the exterior. So the exterior region can effec- 
tively fix the phase in the cooled bubble, suppressing the 
formation of the vacuum states, which would be different 
from that in the bulk liquid. Due to this proximity effect 
the formation of vortices can be also suppressed. 

In both cases of monotonic and nonmonotonic temper- 
ature profile, two mechanisms of the vortex formation are 
important: 

(a) The Kibble-Zurek (KZ) mechanism of the defect 
formation during the quench. For the scenario (ii), where 
the interior region is not separated from the exterior by 
the warmer shell, the KZ mechanism is to be modified to 
include spatial inhomegeneity, which leads to the mov- 
ing transition front. The proximity effect of the exterior 
region is not effective if the phase transition front moves 
sufficiently rapidly 0-^. The modified KZ mechanism is 
not sensitive to the existence of the external counterflow, 
which only role is to extract the formed vortices from the 
bubble. The same KZ mechanism could be responsible 
for the formation of the A-B interfaces, which provides 
another scenario of the B-phase nucleation in the super- 
cooled A-phase 

(b) Instability of the normal-superfluid interface, which 
occurs in the presence of the counterflow Q. 

Here we discuss these two mechanisms (a) and (b) of 
vortex formation during inhomogeneous quench as man- 
ifested in numerical simulations [O. 



II. KZ SCENARIO IN PESENCE OF PLANAR 
FRONT. 

For a rough understanding of the KZ scenario of vortex 
formation let us consider the time-dependent Ginzburg- 
Landau (TDGL) equation for the one-component order 
parameter (OP) ^ = A/Aq: 



To- 



dt 



1 



T(r,t) 



(1) 



Here tq ~ 1/ Aq and are correspondingly the relaxation 
time of the OP and the coherence length far from Tc. 

If the quench occurs homogeneously in the whole space 
r, the temperature depends only on one parameter, the 
quench time tq: 



l-—\Tr 



(2) 



In the presence of a temperature gradient, say, along x, 
a new parameter appears: 



t — x/u 



(3) 



Here u is the velocity of the temperature front which is 
related to the temperature gradient 



UTQ 



(4) 



There exists a characteristic critical velocity Uc of the 
propagating temperature front. At u > Uc the vortices 
are formed, while dit u < Uc the defect formation is either 
strongly suppressed |^,|^ or completely stops [^. 

At slow velocities, u ^ 0^ the order parameter almost 
follows the transition temperature front: 
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T{x - ut) 



, T<Te 



(5) 



In this case the phase coherence is preserved behind the 
transition front and thus no defect formation is possible. 

The extreme case of large velocity of the temperature 
front, ^ oo, corresponds to the homogeneous quench. 
As was found by Kopnin and Thuneberg [|j , if is large 
enough, the phase transition front cannot follow the tem- 
perature front: it lags behind (see Fig. ||). In the space 
between these two boundaries the temperature is already 
below the phase transition temperature, T < Tc, but the 
phase transition did not yet happen, and the OP is still 
not formed, ^ = 0. This situation is unstable towards the 
formation of bubbles of the new phase with ^ 7^ 0. This 
occurs independently in different regions of the space, 
leading to vortex formation according to the KZ mech- 
anism. At a given point of space r the development of 
the instability can be found from the linearized TDGL 
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Order parameter in moving normal/superfluid interface 



from: N.B. Kopnin and E.V. Thuneberg 
Phys. Rev. Lett. 83(1999) 116 
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FIG. 2. The OP distribution from Ref. at nonzero veloc- 
ity u of the planar temperature front in the reference frame of 
the moving front. The higher u the larger is the lag between 
the temperature front at x = 0, where T = Tc, and the OP 
front, which bounds the region with ^ = 0. At > i^c the 
exponentially growing OP fluctuations in the space between 
these two boundaries are not washed out by the moving front 
and lead to the vortex formation according to KZ scenario. 
Dashed line is the Eq.(n[). 



equation, since during the intitial growth of the OP ^ 
the cubic term can be neglected: 



5^ t ^ 



(6) 



This gives an exponentially growing OP, which starts 
from some seed ^auc, caused by fluctuations: 



^(r,t) =^fluc(r)exp 



2rQro 



(7) 



Because of the exponential growth, even if the seed is 
small, the modulus of the OP reaches its equilibrium 
value l^eql = \/l - T/Tc after the Zurek time tz 

tz ^/TQTo . (8) 

This occurs independently in different regions of space 
and thus the phases of the OP in each bubble are not 



correlated. The spatial correlation between the phases 
becomes important at distances where the gradient 
term in Eq. (|^) becomes comparable to the other terms 
at t = tz. Equating the gradient term ^q^^^ ~ (^q/^^)^ 
to, say, the term To^^/Stlt^urek = \/'^o/'?"q^, one obtains 
the characteristic Zurek length scale which determines 
the initial distance between the defects in homogeneous 
quench: 



= Co (rQ/ro)^/^ 



(9) 



We can estimate the lower limit of the characteristic 
value of the fluctuations ^auc = Afluc/Ao, which serve as 
a seed for the vortex formation. If there is no other source 
of fluctuations, caused, say, by external noise, the initial 
seed is provided by thermal fluctuations of the order pa- 
rameter in the volume . The energy of such fluctuation 
is ^^A'^^^Np/Ep, where Ep is the Fermi energy and Np 
the fermionic density of states in the normal Fermi liq- 
uid. Equating this energy to the temperature T ^ Tc one 
obtains the magnitude of the thermal fluctuations of the 
OP 



Ep 



(10) 



Since the fluctuations are initially rather small their 
growth time exceeds the Zurek time by the factor 

Vln|*eq|/|*fluc|. 

The criterium for the defect formation is that the 
time of growth of fluctuations, ~ tz = y'Tgro, is 
shorter than the time t^w = {u) /u in which the transi- 
tion front sweeps the space between the two boundaries. 
Here xo{u) is the lag between the transition temperature 
front and the OP front (see Fig. Thus the equation 
tz = xo{uc)/uc gives an estimate for the critical value 
Uc of the velocity of the temperature front, at which the 
laminar propagation becomes unstable. At large u one 
has xo{u) ~ u^tqTq /4:^Q ^ and thus 



1/4 



^0 V^Q/ 



(11) 



which agrees with estimation i^c = ^v/^z in 0. 

In the case of the neutron bubble the velocity of the 
temperature front is u ~ Ri^/tq^ which makes u ^ 10 
m/s. The critical velocity Uc we can estimate to pos- 
sess the same order of magnitude value. This estimation 
suggests that the thermal gradient should be sufficiently 
steep in the neutron bubble such that defect formation 
can be expected. The further fate of the vortex tangle 
formed under the KZ mechanism is the phase ordering 
process: the intervortex distance continuously increases 
until it reaches the critical size, when the vortex loops 
are expanded by the counterflow. This reproduces the 
most favourable scenario of the vortex formation with 
the cubic law. 
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Rough scenario of vortex formation 
in neutron experiments 



interpretation of results by 
I.S.Aranson, N.B.Kopnin, V.M.Vinokur, 
cond-mat/9905286. 



Superfluid region with Vs=0 is formed 
separated from bulk by vortex sheet 




Two sources of vorticity: 

(1) Vortex sheet 
is unstable to vortex chain. 
(2) Propagating interface 
gives Kibble-Zurek vortices 

Interaction of 2 mechanisms: 

vortex chain screens superflow 

Vs=0 

No Magnus force on 
Kibble-Zurek vortices: 
they decay before the 
vortex chain escapes to bulk 



FIG. 3. Rough scenario of instabiUty of the super- 
fluid-normal interface in the presence of external superflow. 



from: 

I.S.Aranson, N.B.Kopnin, V.M.Vinokur, 
cond-mat/9905286. 




only vortices from vortex layer 
escape into bulk ^He-B 
Number of formed vortices 
-1 

instead of experimental 
N~(f )3-l 

FIG. 4. 3D numerical TDGL simulations of vortex forma- 
tion after the heated region is cooled down (from Ref. 



III. INSTABILITY OF NORMAL/SUPERFLUID 
INTERFACE. 

Another mechanism of the vortex formation has been 
recently found in 3D numerical simulations in Ref. 
It is related to the instability of the normal-superfluid in- 
terface in the presence of the superflow. Let us consider a 
simple hand-waving interpretation of such an instability. 
The process can be roughly splitted into two stages (see 
Fig. I). 

At first stage the heated region of the normal liquid 
surrounded by the superflow undergoes a superfluid tran- 
sition. The transition should occur into the state with 
the lowest energy, which corresponds to the superfluid at 
rest, i.e. with = 0. Thus there appears the superfluid- 
superfluid interface, which separates the state with su- 
perflow (outside) from the state without superflow (in- 
side). Such a superfluid-superfluid interface with tan- 
gential discontinuity of the superfluid velocity represents 
a vortex sheet by definition. Such vortex sheet, on which 
the phase of the OP is not determined, was suggested by 
Landau and Lifshitz for He-II to describe the superfluid 
state of ^He under rotation (see also |l4j and [p!5[). 

The vortex sheet is unstable towards breaking up into 
a chain of quantized vortex lines. The development of 



this instability represents the second stage of the process. 
In numerical simulation the resulting chain of vortices is 
clearly seen (see Figs. ^ and ||). 

The evolution in Fig. |l| is thus caused by the hydro- 
dynamic instability of the normal/superfluid interface in 
the presence of the tangential flow. Since vorticity is 
quantized, such instability leads to the formation of the 
vortex chain only above the treshold required to achieve 
the circulation quantum from the tangential superflow. 
If the counterflow is large the number of vortices in this 
chain TV ^ \vs—Vn\Rh/i^^ i-e. one has a linear law instead 
of cubic. 

Nucleation of the KZ vortices due the motion of the 
superfluid/normal interface is also observed in numerical 
simulations in Ref. [|l^ . It occurs during shrinking of the 
interior region with normal fluid. 

IV. DISCUSSION. 

Two mechanisms of the vortex formation have been 
identified in numerical simulations [Q: (a) vortices are 
formed behind the propagating front due to KZ mecha- 
nism, as discussed in Refs. and in addition (b) vor- 
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Numerical simulation of vortex formation 
in neutron experiment 



I.S. Aranson, N.B. Kopnin, V.M. Vinokur, 
cond-mat/9905286. 




Kibble-Zurek 
vortices 



Kibble-Zurek vortices decay 
because superflow is screened 
by vortex sheet 



Chain of vortices 
(vortex sheet) 



Vortices from vortex sheet 
grow due to Magnus force 
from external superflow 



FIG. 5. 2D simulations, which assume an axisymmetric de- 
velopment (from Ref. [^). 



tices are formed due to the corrugation instability (vortex 
sheet instability) of the front in the presence of external 
superflow. Each of these mechanisms can be derived ei- 
ther analytically for a simple geometry, or understood 
qualitatively with simple physical picture in mind. The 
AKV calculations actually showed that each mechanism 
is fundamental: it does not depend much on the geome- 
try and on parameters of the TDGL equation. Probably 
both mechanisms hold even if TDGL theory cannot be 
applied. 

The interplay of the two mechanisms must depend on 
details of the microscopic physics. In their calculations 
based on TDGL model, AKV found that the chain of 
vortices formed in the process (b) screens the external 
superflow very effectively. The KZ vortices formed in the 
process (a) cannot grow: they decay before the screen- 
ing chain escapes to the bulk liquid. Thus in the AKV 
scenario only the chain of vortices survives. This gives 
the linear dependence of the vortex number N on the 
counterflow Vg — Vn instead of the observed cubic law. 

This does not exclude the possibility of another regime, 
where KZ vortices have enough time to escape to the 
bulk. This is probably what the cubic law found in 
Helsinki experiments tells us. Maybe the latter regime 



cannot be obtained in the TDGL scheme and one must 
discuss the combined dynamics of the OP and quasipar- 
ticles. 

In conclusion, in the period between LT-21 and LT- 
22 the principles of defect formation in inhomogeneous 
phase transition have been developed. 

I am indepbt to N. Kopnin and M. Krusius for dis- 
cussions and I. Aranson and E. Thuneberg who kindly 
provided me with Figures from their papers. 
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